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REDUCTION METHOD WITH SYSTEM ANALYSIS FOR 
MULTIOBJECTIVE OPTIMIZATION-BASED DESIGN 
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Hampton, Virginia 23681 

ABSTRACT 

We present an approach for reducing the number of variables and constraints, which 
is combined with System Analysis Equations (SAE), for multiobjective optimization-based 
design. In order to develop a simplified analysis model, the SAE is computed outside an op- 
timization loop and then approximated for use by an optimizer. Two examples are presented 
to demonstrate the approach. 
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1 Introduction 


There are mainly two classes of methods for handling optimization-based design of complex 
engineering systems. In the context of this paper, a complex system entails two or more 
of the following features: (i) has computationally costly SAE, (ii) has a large number of 
variables and nonlinear constraints, (iii) is multiobjective, and (iv) is decomposable into a 
number of subsystems which hierarchically or nonhierarchically interact with one another. 

The first class of methods are mostly applied to handle features (i) and (ii) (Vanderplaats, 
1984, where further references can be found), and more recently (iii) (see, for example, Zhou 
and Tits, 1993a). In majority of the methods in this class, the system analysis is performed 
outside an optimization loop; often the system analysis (which is a part of an outside loop, 
see Figure 2 for an example) is approximated for use in the optimization loop (the inside 
loop, see Figure 2) in an attempt to reduce the number of costly and detailed analyses. The 
number of variables can also be reduced by design variable linking (Vanderplaats, 1984), this 
is usually done by problem-dependent assumptions (e.g., symmetry in structural design). 
Finally, the number of constraints can be reduced by employing only a significant subset of 
constraints (the active set) in the optimization loop (Zhou and Tits, 1993b; Gill et al., 1982). 

The second class of methods are almost entirely applied to handle single-objective prob- 
lems with the features (i), (ii), and (iv) ( see, for example, Sobieski, 1992; Wu and Azarm, 
1992; Zhao and Azarm, 1993, where further references can be found), with the exception of 
Haimes et al.(1990) method which is applicable to hierarchical multiobjective problems, and 
that of Azarm and Escheuauer (1993) which does not handle the SAE as formulated here. In 
majority of the methods in this class, the system analysis and/or optimization model can be 
decomposed into a number of submodels (subsystems or subproblems). The decomposition 
(which might be performed on the system analysis or optimization model or both) is called 
hierarchic, if it has an overall tree-type structure with two or more levels of subsystems. Ill 
hierarchic decomposition, horizontal interaction in between the subsystems located at the 
same level is not permissible. On the other hand, the decomposition is called uonhierarc.hic, 
if there is no restriction on the interaction which might exist in between the subsystems. 

So far, none of the methods reported in the literature can simultaneously handle all of the 
above-mentioned features which are becoming increasingly common in the design of complex 
engineering systems. In an attempt towards removing this shortcoming, this paper presents 
an extension of a recent work by Azarm and Escheuauer (1993) whereby a multiobjective 
approach for reducing the number of variables and constraints is combined with the SAE 
which might be hierarchically or nonhierarchically coupled. 

The remaining sections of this paper are organized as follows. We present, in Section 2, 
an overview of the definitions and formulation of the problem. This is followed by Section 
3, whereby the solution approach is presented. The method is then demonstrated, as shown 
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in Section 4, via two examples: a simple explosive actuated cylinder (to demonstrate the 
solution steps) and a fairly complex dual-wheel excavator. The paper is concluded in Section 
5 with the final remarks. Parameters and variables for the formulation and examples are 
defined in the nomenclature of Section 8. 

2 Formulation 

As it was stated in Azarm and Eschenauer (1993), the overall multiobjective optimization 
problem is converted into a minmax form in that the collective objective is to minimize the 
maximum loss of all objectives (see also, Osyczka, 1984). Hence, the overall multiobjective 
optimization formulation is given as follows: 

minmax f(z) 

s.t. : g{z) <0 (1) 

where the objective and inequality constraint functions, / and g components, are smooth 
functions of design variables, z. The vector of design variables, z, is partitioned into three 
groups: x, X s , and Y . x represents an ACvector of primary variables, their impact on the 
design optimization is assumed to be global (e.g., variables that contribute to many design 
parts, disciplines, or nonphysical entities). X s represents an .9-vector of secondary variables, 
their impact on the design is assumed to be local (e.g., variables that contribute to detailed 
dimensions or specifications of design parts, disciplines, or nonphysical entities). In general, 
as will be observed in the solution steps (Section 3), the primary variables are included into 
the solution process directly, while the secondary variables are considered indirectly (i.e., via 
step-size variables). Y represents a vector of state variables obtainable as a solution of a set 
of simultaneous (coupled) equations which can be partitioned, for example, into Y , Y , Y , 
such that: 

Y a = f a {x,X„Y b ,Y*) (2) 

Y b = f b (x,X 3 ,Y',Y a ) (3) 

Y‘ = Mx,X 3 ,Y\Y b ) (4) 

The above set of equations represents the SAE. Each equation, also referred to as Contribut- 
ing Analysis (CA), may represent a particular engineering discipline or a distinct physical 
part (a subsystem) of a system, or both (Sobieski, 1990). Note that the division of SAE 
into CAs may: (i) correspond to physical boundaries present in the problem that suggest 
its separation into smaller parts, or (ii) be purely formal as in dividing a set of equations 
into several subsets. These two ways of the SAE partitioning are referred to as physical and 
nonphysical in the paper. 
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As it was stated before, presumably the system analysis is performed outside where the 
objective and constraint functions (optimization formulation) are evaluated. Furthermore, 
it is assumed that in the optimization formulation, the number of variables and constraint 
functions (and perhaps the number of objective functions as well) during the solution pro- 
cess has to be reduced. The reduction in the number of variables is achieved through the 
secondary variables, X s . This has recently been developed for a multiobjective case by 
Azarm and Eschenauer (1993). Briefly, it involves partitioning of secondary variables, X a , 
into several groups: X St , i = 1, • • • , I. Each X Si represents a S\ -vector of secondary variables 
in an optimization subsystem i (which might be an engineering discipline, a physical part, 
or a nonphysical entity). Note that: S = J2i=i $i] the total number of secondary variables 
is equal to the sum of the number of secondary variables in each optimization subsystem: 
i — 1, • • • , /. The vector of secondary variables X Si is replaced by a scalar step-size variable 
Si as follows: 


X 


7iew 

s t 



+ Si di 


( 5 ) 


where di (as given in Section 3) is a descent direction in the subspace of active constraints 
which is obtained from the optimization subsystem post-optimality sensitivity analysis. In 
general, the optimization subsystems (or system) can be different from the analysis subsys- 
tems (here we use the terms: optimization subsystem and analysis subsystem to distinguish 
them from each other). 

In order to reduce the number of constraints, the cumulative function (or the KS function) 
introduced by Kreisselmeier and Steinhauser (1979), which was also used by Sobieski et al. 
(1985), is utilized here. As an example, g^s represents a cumulative function of <7i, • • • , gj: 


9h 


= -In 


J2 ex P(Pl h) 

U=i 


( 6 ) 


where p is a user controlled coefficient. Sobieski et al.(1988) were first to develop a technique, 
based on the KS function, to convert a multiobjective optimization problem into a single 
objective form. The KS function can also be applied to convert a multiobjective optimization 
problem formulated in a minmax form into a single objective form. Figure 1 shows an 
example of this whereby /* s and g^ s are the KS functions for the objectives and constraints, 
respectively. Note that, as shown in the Figure 1, depending on the value of p there 
might be a gap in between the minimum of fk s and the minmax(fi, / 2 , fa), points B and 
A, respectively. These two points however will coincide when p — ► oo. Also, note that the 
choice of the KS function as a means to reduce the problem to a single objective form implies 
that the objective that has the steepest slope after normalization controls the optimization. 


Finally, the SAE can be approximated in the optimization formulation in which the 
number of variables and constraints are reduced as well, i.e., the optimization problem of (l) 
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Figure 1: KS Function for Objectives and Constraints 

can be rewritten as follows (x and $ are varied while Y fixed): 

minmax f (x,.Si, ■ • • , .s/, F) 

s.t. : gks (®, «!,•••, 5;,?) < 0 (7) 

where Y € [F[, Y represents a linear approximation of the Y which is computed by 
the SAE, Y° t and Y° u represent lower and upper move limits on the current value of Y. The 
advantage of (7), when compared with (1), is that the number of constraints is reduced to 
one and the number of variables is reduced from (N+S'+number of state variables) in (1) to 
(Af -|- /) in (7). As an example, in a large-scale structural optimization problem, the reduction 
in the number of variables and constraints can (among other reasons) substantially reduce 
the computer storage for an otherwise very large jacobian matrix needed by the optimizer. 

3 Solution Steps 

The solution steps are summarized as follows (see also the flow-chart in Figure 2): 

• Step (0): Identify primary, secondary, step-size and state variables. Initialize with 

x = x°, X 3l = , Si = 0 and set k=0, 

• Step (1): Solve system analysis equations for V , and perform system sensitivity analysis 
(GSE as described by Sobieski, 1990) to obtain f) Y /dx and /0X S . 
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Figure 2: Flow-Chart of the Solution Steps 


Step (2): if k=0 then solve (8) where x is varied while Y is fixed: 

minmax f(x, Y) 

*.t. : gks(x,Y) < 0 

otherwise (i.e., if k > 1), solve (9) where x and s are varied while Y fixed: 

minmax f(x, .si, • • • , .s/, F) 
s.t.: 9ks(x,Si,- ■ • < 0 

where in (8) and (9), there are move limits on Y , i.e., Y € [ Y h Y u ]. 

Step (3): As a part of post optimality sensitivity analysis of (8) or (9), compute: 







and find: 


X"' w = X°[ d + .s, d, (11) 

• Step (4): Set k=k+l, repeat steps (l)-(3) until convergence is achieved. 

In the above-mentioned solution steps, u and A are Lagrange multipliers corresponding to 
the objective and constraint functions of (8) or (9). Furthermore, Y is obtained as follows. 

? = Y + ^{{dY Idxfi&Xi) + £{(< dY/dX Sl )(dX Sl /dsi)Asi } (12) 

t * 

— C- ~ 

where Y is the current value of Y . Also: 

(d{f or g ks )/dx) = (d(f or fii*.)/^) loca l + E W or 9ks)/dY { )(dY 702) (13) 

t=a,6,c 

(d(f or gks)/d$i) = ^ {(dY* /dX s )(d(f or gks)/dY')}(dX s /dsi) (14) 

t'=a,6,c 

The minimax problems of (8) and (9) are solved by the subroutine FSQP (Zhou and Tits, 
1991). FSQP (Feasible Sequential Quadratic Programming) is a set of Fortran subroutines 
which implements algorithms which are described and analyzed by Panier and Tits (1993), 
Bonnans et al. (1992), and Zhou and Tits (1993a). 


4 Examples 

Two examples are presented here to demonstrate the method developed in the paper. Both 
examples are selected from the literature where they are formulated as a single-objective 
problem. They are revised here to form multiobjective examples. 

4.1 Example 1: Explosive Actuated Cylinder 

This example is constructed from a well-known single-objective optimization problem, a 
minimum length design (/i) of an explosive actuated cylinder (Papalambros and Wilde, 
1979). The constraints for this example express the specifications for: kinetic energy (</i), 
wall stress (g 2 ), and geometry constraints {gz-gs). Its formulation is revised here by taking 
out one of its constraints (piston force) in the original formulation and including it as an 
additional objective (f- 2 ) in the problem. As shown below, both objectives have been scaled 
so that they are of the same order of magnitude: 

minmax { /i , /a } 

/l = { Z 1 T Z 2 )f L in ax 
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i 
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i 1 mi "i !»t'l Mini H 'i 


f‘i = [( 1000 ^/ 4 )^ 4 ^ 5 ]/ F max 

s.t. : (15) 

<7i : (Wm*«(l - 7))/[1000 z 4 i> 7(4 1_7) - u i (1_7) )] - 1 < 0 

fi = v c + {n/4)zizl 
v-i = Vi + (7t/4)z 2 ^5 

<72 : a e N/a Y - 1 < 0 

&e = (cr'i - cr x a 2 + 0 2 ) 1/2 

<Ti = * 4(4 + *5 ) / ( *3 - *5) 

<r 2 = *4 


<73 : 

2.3/ Dmax - 1 < 0 


<74 : 

(*1 H" *2)/ L m(l x 

1 < 0 

<7s : 

*5 < *3 



V 

0 

. 

II 

•,5 


For convenience, the nomenclature for this example is given in the appendix (see also Pa- 
palambros and Wilde, 1979, for further details). 


4.1.1 Solution 

Following the solution steps in Section 3, we assume: x = (-* 1 , 24 )', X 3 - ( 22 , 23 , 25 )', Y = 
(Y a ,Y b Y = (v u v 2 y. We then initialize with x° = (0.4,0.34)', X° = (1.2,0.82,0.67)', .s°=0, 
and k=0. Next, as in step ( 1 ) of the solution steps, we select and partition the SAE into 
two nonphysical analysis subsystems: subsystem a and subsystem b (Figure 3 shows the 
interaction between these two subsystems). We then form the global sensitivity equations 


a 2 

Y = v = v +TT/4 z „ z . 

1 c 15 

subsystem a subsystem b 



b 2 

Y = v = v +7174 z z 


2 1 2 5 


Figure 3: Analysis Subsystems for Example 1 


(Sobieski, 1990): 

( 1 —dv x /dv 2 \ / dvi /dz k ) _ f df a jdz k \ 

\ -dv-i/dv 1 1 j ^ dv 2 /dz k ) \ df b /dz k ) 

where: 

f a = U c + (?r/4)2 1 25 


A: = 1,2,5 (16) 

(17) 
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Test Cases 

a 

b 

c 

d 

e 


HOT 

0.722 

0.723 

0.724 

0.724 


HOT 

0.722 

0.723 

0.724 

0.724 


HD 

H 

H 

HOI 

0 


wm : $ 


na 

mi 

1 .45 


Bpii 


0.97 

msEi 

0.89 

za 

23.4 

23.4 

23.04 

23.04 

22.5 

*5 

0.166 

0.166 

0.167 

0.167 

0.169 


Key: 

a: original double-objective problem 
b: original double-objective problem was 

reduced to single-objective by KS function 
c: original multi-constraint problem was 

reduced to single-constraint by KS function 
d: no. of variables and constraints are reduced 
e: no. of variables and constraints are reduced 
with linearization on <71 

Table 1: Summary of Results for Example 1 


f b = uj + (ir/A)z 2 zl (18) 

As an example, OSE for 25 will be as follows: 

+ 1 0 \ ( dvxjdzs \_( {it 1 2) z\Zs \ nq) 

-1 1 ){dv 2 /dz 5 )-{(n/2)z 2 z 5 ) 1 ^ 

With k=0 at step (2), we formulate (8) with p=150, and solve for updated primary variables: 
£=(0.75,2.3)*. We also obtain (as a by-product of FSQP): (u>i,u> 2 )= (0,1) and A=0.2. Next, 
we go to step (3), whereby (10) we compute d u d i) d 3 for X], X s 2 , X ,, respectively. We then 
express the secondary variables via (11) as a function of the step-size variable $ (new values 
for the step-size .s and primary variables x are then obtained via (9) for k > 1). We now 
go to step (4) and set k=l and repeat steps (l)-(3). This iterative process continues until 
there is no further improvement in the value of objective functions. A summary of results 
are given in Table 1. As shown in the table, the results from the solution steps (test case e) 
compare well with those obtained when the original problem is solved (test case a), or when 
the number of objectives are reduced to one (test case b), or w'hen the number of constraints 
are reduced to one (test case c), or when the number of variables ( including the step-size 
variable s) are reduced to three and the number of constraints to one (test case d). 
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4.2 Example 2: Dual- Wheel Excavator 

This example (Figures 4 and 5) is constructed from a single-objective optimization problem 
(Wilson, 1992). It involves a minimum weight design of a 120-inch threaded hub-and-shaft 
assembly for a dual-wheel excavator. There are 9 variables and 25 constraints in this example. 
The constraints in the example express specifications on: hub stress (<?i, ( 72 ), stress on shaft at 
various critical positions (< 73 -< 7 is), and other practical constraints ( gwg-is )• Here, the original 
single-objective problem (sum of the weights for the hubs and shaft) has been converted into 
two objectives: weight of the shaft (/ 1 ), weight of the hub (/ 2 ), with both objectives scaled 
to be of the same order of magnitude: 


minmax {/i,/ 2 } 

f\ = {psff/4){Zspz 3 -f Zsh z 4 + 2 Zb D% + 2z&z\ 

+2 z 7 z] + 2 Z c z\ + 2z 9 «|}/3300 

f-2 = (psft/ 'ty{[ z r{( z \ + «) 2 — z \)] + [("7 + 5 2 )/2(Z)j 1 

-(*1 + «)' 2 )] + K*i + B,)[2{D\ 0 - D 2 u )}/mo 

s.t. : 


9 \{ z \i 

^7) 


Gahj Gy "h &mhj Gy 

//. < 0 

<72(21, 

27) 


&ah/&e &mh/ G u t 

1//. < 0 

<73(23, 

27,2 

s) : 

Gab/ Gy “H &mbl Gy 1 / fs ^ 0 

<74(23, 

27,2 

s) : 

Gab/ G e “ 1 “ &mbj &ut 

l //.<0 

<75(21, 

*t) 


Gauf Gy + @muj Gy 

L //.<0 

< 7 e (21, 

*t) 


Vau/Ve + &mul&ut 

l// s < 0 

<77(21, 

*r) 


G athd j Gy H" &mthd/ Gy 

- 1 //, < 0 

< 7 s( 2 l, 

27): 


Gathdj &e H~ G 'mthdf Gut 

- 1 //, < 0 

<79(21) 

* 


G axthj G y + & mxth! Gy 

- i/f. < 0 

<710(21 ) : 


G axthj G e + &mxthj G ut 

- 1 If. < 0 

< 7 n( 2 i, 

- 5 ? * 

-<0 : 

& (tend/ &y H” &mendj&y 

- 1 /fs < 0 

<712(21, 

- 5 , * 

*): 

Gaend j G e ” 1 " &mendj G u 

- 1 //, < 0 

<71.3(21, 

-’5) 


G dent/ Gy “h ^7? lent j Gy 

- 1 ffs < 0 

<714(21, 

-5) 


Gaent / G e G men tj Gut 

- 1 Ifs < 0 

<715(21, 

-2) 


Gare vj Gy “ 1 “ G inT pyj Gy 

- 1 If, < 0 

< 7 ie( 2 i, 

*2) 


Garev j G e ”h G mrev j<J u i 

-l //.<0 

<717(21, 

25, 2 

<9): 

Gats/ Gy 4 “ G ines j Gy 

1 If, < 0 

< 7 is( 2 i, 

25,2 

■9) : 

Gats/Gf. + G mes /G u i — 

• Ilfs <0 


( 20 ) 
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<t» z 6 <t> z 5 4> z l 0 Z 2 < t )D B 4> z 4 4> z 3 



Figure 5: Shaft of the Dual-Wheel Excavator 


g™(zs) 


0.127 N s - 2 8 < 0 

<720 (Z2, 


z 2 - 0.95 D b < 0 

<721 (z 2 , 

- 3 ) 

0 . 8 D# — z 2 < 0 

< 722 ( 2 . 3 , 

24 ) 

Skidr + Dfl /2 — 3 : 4/2 < 0 

< 723 ( 25 , 

2 6 ) 

*» + 0.0165 -* fl <0 

< 724 ( 21 , 

Z 'i ) 

Zi + 0.0267 - z 2 < 0 

< 725 ( 21 , 

*.) 

0.25z, - z 5 < 0 


For convenience, the nomenclature for this example is given in the appendix (see also Wilson, 
1992, for further details). 

4.2.1 Solution 

Following the solution steps in Section 3, we assume: x=(zi, z 7 )', X Si =(z 2 , 2 3 , 24 , zs) 1 , 
X, 2 = ( - 5 , ^ 6 ) 2 g)h Note that the secondary variables are grouped (optimization subsystems) 
according to the physical makeup of the excavator; X Sl represents shaft dimensions in be- 
tween the two hubs while X S2 represents those outside the two hubs (Figures 4 and 5). 
State variables Y, or the SAE (Figure 6 ), for this example are established so that they only 
contribute to the hub stress constraints (<71 and g 2 ). The hub is modeled (Wilson, 1992) as a 



Figure 6: Analysis Subsystems for Example 2 
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ring plate of linearly varying thickness subject to concentrated transverse load and bending 
moment at its outer edge. The SAE is partitioned into three nonphysical analysis subsys- 
tems: subsystem a ($'“), subsystem b (V" 6 ), and subsystem c. (f c ). Note that, the analysis 
subsystems as shown in the Figure 6 are different from the optimization subsystems. The 
SAE for subsystem a represents: 


subsystem a: 


V? = n = /T(*i,*t) 

Y b = A = /a (n) 


( 21 ) 


where n is the hub’s thickness constant, and A is an eigenvalue computed iteratively (as shown 
by an arched arrow on subsystem a in Figure 6) by the following characteristic equation. 
The characteristic equation is obtained from partial differential equations established for 
the hub based on the excavation loading conditions and using small-deflection plate theory 
(Conway, 1958; Wilson, 1992): 


A 4 - 2A 3 (n + 2) + A 2 [(n + 2) 2 + n(l - v) - 2m 2 ] 

+ A(n + 2) [2m 2 - n(l - i ')] + m 2 [m 2 - 1 - (3 + un){n + 1)] = 0 


( 22 ) 


subsystem b: < 


A is obtained for each to (only 11 terms are considered here, i.e., m=0 to m=10) of the de- 
flection {w) equation in subsystem b. As shown below, subsystem b computes the deflection 
and loads (shear force, Q r , and moments, M rt M e , M rS ) on the hub: 

Y b = W = T,m=o( Ar * 1 + + ClM + Dr **) cos ( m6 ) 

Y 2 b = Qr = f${w,z u z 7 ,n) 

Y b = Mr = f b i{w , Z\ , 27 , n) 

Y 4 b = M 8 = f%(w,z u z 7 ,n) 

Yf = M r$ = f 4 {w,z u z 7 ,n) 

Where: Q r = D(r)[d/dr}((d 2 w/dv 2 ) + (1 /r)(dw/dr) + (1 /r 2 ){d 2 w/d9 2 )) 

M r = —D(r)[(d 2 w/dr 2 ) + v((\/r)(dw/dr)+(l/r 2 )(d 2 w/d0' 2 ))] (24) 

Me = -D{r)[{l/r)(dw/dr) + (l/r 2 )(d 2 w/0 2 ) + v(d' 2 w/dr 2 )} 

M r o = -(1 - u)D(r)[{\/r){d 2 w/drdO)] 

and D(r), the flexural rigidity, is a function of hub radius, r, and 8 is an angular coordinate 
on the hub. 

Finally, subsystem c. represents the boundary conditions which can be used to find 
A,B,C,D coefficients for the above-mentioned deflection ( w ), Wilson (1992): 


subsystem c: 


Y{ = A; Y 2 c = B\ Yf = C\ Y 4 = D 
w = /{(A, B,C, D, A) = 0 

dwfdr =/l(A,B,C,/9,A)=0 

^(forces) = /|(A, C, D , A, n, Q r , M r e) = 0 
^(moments) = /J(A, B, C, D, A, n, M r ) = 0 


(25) 


Initially, we assumed: x=(12.1, 8.9) ! , X S] — (14.4, 16.6, 19.1)‘, X„ 2 (9.2, 12, 1.8) , .Si 0, 
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Test Cases 

a 

■1 

■B 

US 

— n 

/1 

0.80 

0.79 

fm 

Hal 

l 

h 

0.80 

0.80 

0.80 

0.81 


z\ 

1EE1 

12.7 

■rail 

11.9 

12 

Z'2 

log 


ca 

13.6 

13.9 

-3 



lUl 

15 

14.8 

z 4 

17.6 

nm 

Bl 

18.6 

18.2 

Z 5 


Bltl 



9 

z 6 

■9 


fern 

11.8 


zr 

■29 


5.9 

6 

H| 

Z S 

9B9 

9B9 


12 


Z 9 

1.3 


1.4 

1.2 

1.2 


Key: 

a: original double-objective problem 
b: original double-objective problem was 

reduced to single-objective by KS function 
c: original multi-constraint problem was 

reduced to single-constraint by KS function 
d: no. of variables and constraints are reduced 
e: no. of variables and constraints are reduced 
with linearization on <71 and <72 

Table 2: Summary of Results for Example 2 

•S2 == 0> and k=0. We then followed the solution steps in Section 3 to obtain the results which 
are summarized in Table 2. Finally, it should be stated that the minimax solution reported 
here (for case a of Table 2, the total hub-and-shaft assembly weight is 10,640 lb) is different 
from that reported by Wilson (1992) in which the problem is formulated in a single objective 
form (sum of the weights for the hubs and shaft was 10,211 lb). These two solutions are 
essentially two different Pareto solutions for the example as formulated here (Osyczka, 1984). 


5 Concluding Remarks 

In this paper we have discussed a method for combining a reduction technique with system 
analysis equations for multiobjective optimization problems. The main characteristics of the 
method presented are that: (i) the number of variables and constraints can be reduced, (ii) 
the SAE, which might be costly and hierarchically or nonhierarchically decomposed, is per- 
formed outside the optimization loop, and (iii) it is applicable to multiobjective optimization 
problems. 

The method has been demonstrated by two examples: (i) a simple explosive actuated 
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cylinder, and (ii) a fairly complex dual-wheel excavator. It has been shown that, foi both 
examples, when the reduction measures (test cases (b)-(e) for both examples) are applied 
they can usually obtain a solution fairly close to that of the original problem (test case 
(a)). The small difference in the solution could be eliminated, for example, by increasing the 
number of variables, constraints, etc., to resolve a lack of sufficient degrees-of-frcedom. 
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8 Nomenclature 


Formulation: 

(CL4) rt 

d x 

l 

fks 

9 

9ks 

I 

k 

N 


S 

S ' , 

x 

x s 

Xs, 

Y 

ya 

- c ~ c 

Y„Y« 

Y 


= Contributing analysis equations in analysis subsystem a 
= Descent direction in optimization subsystem i 
= Vector of objective functions 
= Cumulative objective function 
= Vector of inequality constraint functions 
= Cumulative constraint function 
= Number of Optimization subsystems 
= Iteration counter 
= Number of elements in tildex 
= Step-size variable in optimization subsystem i 
— Number of elements in X s 
= Number of elements in X Si 
= Vector of primary variables 
= Vector of secondary variables 

= Vector of secondary variables in optimization subsystem i 
= Vector of state variables 

= Vector of state variables in analysis subsystem a 
= Lower and upper move limits on the current value of Y 
= Linear approximation to Y 
= Vector of overall design variables 
= Lagrange multipliers for constraint functions 


15 



max 

max 

max 


p 

U) 

Example 1: 

D 
F, 

L 
N 

v u v 2 

Vc 

Wmin 
S\ 

~2 
~3 
Z 4 
~5 

7 

0Y 

01, <*2 


= User defined coefficient for KS function 
= Lagrange multipliers for objective functions 

= Maximum allowable cylinder outside diameter, 
= Maximum piston force, 700 lb 
= Maximum cylinder total length, 2 in 
= Safety factor, 3 

= Initial, final volume of combustion, in 3 
= Fixed chamber volume, 0.084 in 3 
= Minimum kinetic energy, 600 lb-in 
= Unswept cylinder length, in 
= Working stroke of piston, in 
= Outside diameter of cylinder, in 
= Initial pressure of combustion, ksi 
= Piston diameter, in 
as Ratio of specific heats, 1.2 
= Equivalent stress, ksi 
= Yield strength, 125 ksi 
= Principal stress, ksi 


in 


Example 2; 

a 

B\ , B 2 
D\o 
On 
D b 

fs 

N s 

Shldr 

*f - 
**\ 

Z B 

Z c 

Z$H 

Zsp 

A 

V 

PS 

&ab 
& aend 
&aent 
&ae$ 

&ak 

&athd 

&au 


= Diameter constant for the hub, in 
= Hub dimensions, in 
= Hub dimension, in 
= Hub dimension, in 
= See Figure 5, in 
= Safety factor 
= Number of seals 

= Radial height of the shoulder at the bearings, in 
= Shaft dimension {i = 1, - ■ • , 9), see Figure 5, in 
= See Figure 5, in 
= See Figure 5, in 
= See Figure 5, in 
= See Figure 5, in 
= Eigenvalue 
= Poisson’s ratio 
= Density, lbm/in 3 

= Alternating stress on shaft at bearings, ksi 

= Alternating compressive stress on shaft at end connections, ksi 
= Alternating tensile stress on shaft at end connections, ksi 
= Alternating shear stress on shaft at end connections, ksi 
= Alternating stress on hub, ksi 

= Alternating stress on shaft due to reverse loading, ksi 
= Alternating stress on ACME threads, ksi 
= Alternating stress on shaft at ACME thread undercut, ksi 


% 
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Caxth = Alternating shear stress on ACME threads, ksi 

cr e = Endurance limit, ksi 

crjnb = Mean stress on shaft at bearings, ksi 

&mend = Mean compressive stress on shaft at end connections, ksi 

Cment = Mean tensile stress on shaft at end connections, ksi 

cr m es — Mean shear stress on shaft at end connections, ksi 

Vmii = Mean stress on hub, ksi 

Vmrev — Mean stress on shaft due to reverse loading, ksi 

&mthd = Mean stress on ACME threads, ksi 

v mu = Mean stress on shaft at ACME thread undercut, ksi 

v mxth = Mean shear stress on ACME threads, ksi 

cr u t = Ultimate tensile strength, ksi 

v y = Yield strength, ksi 
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